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Schrodinger models for solutions of the Bethe-Salpeter equation in Minkowski space 
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By application of the 'geometric spectral inversion' technique, which we have recently generalized to 
accommodate also singular interaction potentials, we construct from spectral data emerging from the 
solution of the Minkowski-space formulation of the homogeneous Bethe-Salpeter equation describing 
bound states of two spinless particles a Schrodinger approach to such states in terms of nonrelativistic 
potential models. This spectrally equivalent modeling of bound states yields their qualitative features 
(masses, form factors, etc.) without having to deal with the more involved Bethe-Salpeter formalism. 

PACS numbers: ll.10.St, 03.65.Pm, 03.65. Ge. 

Keywords: Bethe-Salpeter equation, nonrelativistic potential model, geometric spectral inversion 



!>• . I. INTRODUCTION: MOTIVATION AND INCENTIVE 



Within the framework of relativistic quantum field theory, the appropriate tool for the description of bound states is, 
Q^j in principle, the Bethe-Salpeter (BS) formalism PHI]- In this approach, a bound state B(P), of momentum P and mass 
M, is described by its BS amplitude, which in configuration-space representation is defined by the matrix element of the 
time-ordered product of the field operators of all bound-state constituents between vacuum |0) and bound state |B(P)) . 
In momentum-space representation, the BS amplitude, upon splitting off the center-of-momentum motion of the bound 
state and suppression of all indices generically denoted by <&(p, P), encodes the distribution of the relative momenta p 
of the bound-state constituents. It satisfies a formally exact BS equation involving two kinds of dynamical ingredients, 
namely, for bound states composed of n constituents, (a) the propagators Si(pi) (i = 1, 2, . . . , n) of the n constituents of 
\ respective individual momenta pi and (b) its BS interaction kernel K, a fully truncated (2 n)-point Green function of the 
f — ■ n bound-state constituents, perturbatively defined as the sum (of the countable infinity) of all 'BS-irreducible' Feynman 
graphs for n-particle into n-particle scattering. For two bound-state constituents, the BS equation is of the generic form 

<P(p,P) = j^S 1 (p 1 ) J d i q K(p, q ,P)<S>(q,P)S 2 (-p 2 ) . (1) 

o ' 

Physical considerations provide a profound motivation to formulate one's BS framework in Minkowski space, with the 
t-H I pseudo-Euclidean space-time metric tensor g^ v — diag(+l, —1,-1,-1). In the Minkowski-space formulation, however, 
finding solutions to the BS equation may be heavily impeded by the presence of singularities induced by the propagators 
of the bound-state constituents or its BS interaction kernel. As a remedy, by assuming that analytic continuation of the 
Minkowski-space formalism is possible and the Cauchy integral theorem is applicable, it has been proposed to study the 
BS equation in Euclidean space, with metric = <5 M „, reached by a procedure misleadingly labeled Wick 'rotation' 
The Euclidean-space formulation facilitates making contact with lattice field theory, usually defined in Euclidean space. 
Solutions of the BS equation provide the set of mass eigenvalues M and associated BS amplitudes $ of the bound states. 
The mass eigenvalues arising in Minkowski-space and Euclidean-space formulations of a given BS equation are identical. 
For BS amplitudes, however, complicated analyticity structures of the BS equation in the complex plane cause troubles: 
The BS solutions derived in one formulation may differ from those obtained in the other one. That is to say, the analytic 
continuation to Minkowski space of some solution to the BS equation in Euclidean space may bear no resemblance to its 
counterpart found as solution to the same BS equation in Minkowski space. Since the BS amplitudes determine physical 
observables such as decay constants and form factors, knowledge of the Minkowski-space amplitudes is highly desirable. 

This dilemma between, on the one hand, the comparative ease of deriving solutions to a Euclidean-space BS equation 
and, on the other hand, the need of physical applications for BS amplitudes constituting solutions to a Minkowski-space 
BS equation can be tentatively resolved by developing — of course, only approximately equivalent — Schrodinger models. 
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This may be effected by fixing the interaction potential entering in the Schrodinger Hamiltonian by spectral inversion of 
the bound-state mass eigenvalues M arising from the easier-to-accomplish solution of the Euclidean-space BS equation. 
The wave functions obtained as solutions of the resulting Schrodinger equation allow us to compute decay constants and 
form factors of bound states, provided we succeed in acquiring control of the uncertainties introduced by such modeling. 
Obviously, we may estimate the accuracy of the envisaged Schrodinger models by applying our proposal for resolution of 
the dilemma to some known solutions of the Minkowski-space BS equation in order to extract an associated Schrodinger 
potential and comparing the outcome of the corresponding Schrodinger equation with the findings of the BS framework. 

Recently, renewed attempts of solving the BS equation for two-particle bound states in Minkowski-space formulation 
have been undertaken [5l4l3l|. The basic idea advocated for in Ref. @ is to remove the singularities of the BS amplitudes 
by considering an equivalent integral transform of the BS equation (fT]) obtained by projection onto the light-front plane, 
and to take advantage of a particular integral representation of the BS amplitude P) proposed by Nakanishi 14] , in 
order to obtain a nonsingular integral equation that is straightforward to solve numerically. This route has been app lied 
to bound states consisting either of two identical spin-0 bosons [5l-flC)j or of a spin-i fermion and its antiparticle [ilTfiljl ]. 
Any such bound state arises, in the case of two scalar constituents, from their couplings to a further scalar boson or, in 
the case of fermionic constituents, from their couplings to a scalar, a pseudoscalar, or a massless vector boson; for bound 
states of scalars, the interactions are taken into account in the BS equation by considering the BS kernel either in ladder 
approximation, which amounts to the iteration of single-boson exchanges, or in ladder-plus-cross-ladder approximation. 
Moreover, for the bound-state constituents the studies (5l4l3l| employ, for simplicity, the free-propagator approximation. 

In order to get an idea of the behaviour of such Schrodinger interaction potentials to be expected to arise in the course 
of spectral inversion, we recall that there is a well-paved path of simplifications leading from the relativistic BS equation 
to its nonrelativistic or 'static' Schrodinger reduction. The sequence of necessary steps involves several well-defined and 
thoroughly studied approximations to the BS formalism (for brief reviews of this reduction, consult, e.g., Refs. [IB4I3) : 

1. In some instantaneous limit, realizable if in the bound-state's center-of-momentum frame fixed by P = (M, 0) the 
BS kernel takes the form K(p, q, P) = K(p, q), the BS equation may be reduced to the instantaneous BS equation 
(for attempts in these directions, consult, for instance, Ref. [lSj and references therein) for the Salpeter amplitude 



2. The additional assumption of free propagation of all bound-state constituents with effective masses encompassing 
the dynamical self-energy effects leads to the Salpeter equation jl~Qj ] - (Note, however, that in quantum field theory 
the Dyson-Schwinger equations relate every n-point Green function to at least one (m > n)-point Green function. 
This means, in particular, that the propagators, i.e., the 2-point Green functions, and the n-point Green functions 
entering in the BS kernel cannot be chosen independently: the use of free propagators might be incompatible with 
the feature of confinement exhibited by quantum chromodynamics, the theory describing the strong interactions.) 

3. Dropping all negative-energy contributions simplifies Salpeter's equation to the reduced Salpeter equation [2"oU24j . 

4. Furthermore, ignoring all spin degrees of freedom of all bound-state constituents and assuming the BS interaction 
kernel K to be of convolution type, i.e., to depend only on the difference of the involved relative momenta p and q, 
K(p, q) = K(p—q), yields the spinless Salpeter equation. Therein, the interactions manifest in form of a potential 
arising, in configuration space, as the Fourier transform of this kernel K{p—q). This bound-state equation may be 
viewed as a generalization of the Schrodinger equation towards relativistic kinematics. Concise reviews of various 
aspects and facets of semirelativistic approaches to the bound-state problem may be found in, e.g., Refs. [25U271 ]. 1 

5. In an ultimate static limit, replacing in the latter equation of motion the relativistic form of all one-particle kinetic 
energies by the corresponding nonrelativistic approximation, we eventually end up with the Schrodinger equation. 

Let us begin our analysis by inspecting the simplest case: bound states of two scalar constituents In the ladder 

approximation, the only contribution to the BS interaction kernel derives from single-particle exchange. Apart from the 
couplings of the exchanged particle to the bound-state constituents, the BS kernel is then nothing but the propagator of 
the exchanged particle. For a scalar boson with mass fi, its free propagator is given by S(k) = i (fc 2 — /j 2 ) -1 . The Fourier 
transform in three dimensions of the instantaneous approximation to this propagator, i.e., of its remnant i (fc 2 +/i 2 ) -1 , is 
proportional to the configuration-space Yukawa potential V(r) — — exp(— fj/r) /r, thus singular at the origin r = \x\ = 0. 
As a consequence, depending, clearly, on the proximity of the system described by the BS equation to the nonrelativistic 
Schrodinger limit, the outcome of any spectral inversion may be potentials resembling, to some extent, the Yukawa type 
but modified, of course, by the various effects ignored on the way down to the static limit, such as relativistic kinematics 




1 A related approach is the quasipotential formalism devised by Todorov |2g( . 
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or higher-order contributions to the BS interaction kernel; cross-ladder terms are but the simplest example of the latter. 
Accordingly, we have to devise and utilize an inversion technique that is capable of dealing also with singular potentials. 

To this end, we recently generalized the earlier geometric spectral inversion [20M33I ] to treat singular potentials Js3| ■ 
We suppose that f(r) is the shape of the potential and v > is the coupling parameter in the Schrodinger Hamiltonian 
H = — A/(2m) +v / (r). In this inversion technique, a functional sequence is built which starts from a seed /^(r) and 
reconstructs the potential shape f(r) from a given spectral function E = F(y) that defines how a discrete eigenvalue E 
of H depends on the parameter v. The most relevant earlier paper is Ref. 34], which also includes a proof of uniqueness 
for the inverse for a large class of singular potentials. Here, the inversion sequence is defined in Sect. Ill CI a statement 
of the uniqueness theorem may be found in Sect. Ill Dt The principal goal of the present paper is to take for F(v) the 
solutions to the Minkowski-space BS equation (114 1 Oil 13} for bound states of two scalar constituents, and to reconstruct 
directly from this set of data an effective potential shape f(r) in the Schrodinger model defined by the Hamiltonian H. 

The outline of this paper is as follows. In Sect. [Ill we summarize enough of the geometric spectral inversion theory to 
make this paper essentially self-contained. In Sect. IIIH we present the spectral data F{v) from the BS solutions [a-tlOj. 
In Sect. HYl we apply the functional inversion sequence of Sect. Ill Cl to these data to construct the effective potential 
shape. In Appendix |A| we sketch the nonrelativistic reduction of the Bethe-Salpeter formalism for scalar bound-state 
constituents along a route which mimics to the utmost possible extent the case of fermionic bound-state constituents. 



II. GEOMETRIC SPECTRAL INVERSION 



We consider the discrete spectrum of a Schrodinger Hamiltonian operator 

H = -A + vf(r), r=\\x\\, (2) 

where f(r) is the shape of an attractive central potential, and v > is a coupling parameter. We shall assume that the 
potentials are monotone non-decreasing and no more singular than the Coulomb potential f(r) — —1/r. The arguments 
we use apply generally to the problem in d > 1 spatial dimensions, but, for definiteness, we shall usually assume that 
d = 3. The operator inequality (35l . l36j 



fd/2 - 1 



2 



-A> ) > d>3, (3) 

implies that a discrete spectrum exists for sufficiently large coupling v > 0. For d — 3, the Hamiltonian H is bounded 
below by 



E > min 

r>0 



1 

4r 



■2 ■ v f( r ) 



(4) 



and a simple trial function can be used to establish an upper bound to E. Thus, we may assume, in particular, that the 
ground-state energy may be written as a function E = F(v). An explicit example of the class of problems we consider is 
provided by the Hulthen potential, whose shape is given by f(r) = — l/(e r — 1) and whose s-state (1 — 0) eigenvalues E n 
are given [371 ] exactly for d — 3 by the formula 



2 \ 2 

v — n 



E " = F ^ = -{^~) > v>n 2 , 71 = 1,2,3,.... (5) 

The problem discussed in the present paper may be stated as follows: given, for example, the curve Fi(v), can we use 
this spectral data to reconstruct the potential shape f(r)7 We call this reconstruction a 'geometric spectral inversion'. 



A. Exact representation of spectral functions by kinetic potentials 

The discrete spectra of operators bounded from below, such as H — — A + v f(r), may be characterized variationally 
[38} . Thus, the ground-state energy may be written 

F(v) = inf Ub, Hi/}) . (6) 
11-011=1 

Since H depends on the coupling v, so therefore does the domain T>(H). However, for the problems considered, either H 
has discrete eigenvalues, perhaps for v greater than some critical coupling v\, or the entire spectrum of H is discrete for 
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v > 0. The kinetic potential f(s) associated with a given potential shape f(r) is defined (for the ground state V>) by a 
constrained minimization in which the mean kinetic energy s = (—A) is kept constant: 

f(s)= inf • (7) 

11^11=1 

The eigenvalue F(v) of H is then recovered from f(s) by a final minimization over s: 

F(v)=min[s + vf(s)]. (8) 

The spectral function F(v) is concave (F"(v) < 0); moreover, it has been shown [29] that 

F"(v)f" (,)= * <o. (9) 

Hence, -F(w) and /(s) have opposite convexities and are related by the following Legendre transformations / O F [39| : 

/(«)=*"(«), a = F(i;)-i;F'(v) , (10) 

- = -/», ^ = /»-*/»• (11) 

_F(u) is not necessarily monotone, but the kinetic potential f(s) is monotone decreasing. Equation (|10p enables us also 
to use the coupling as a minimization parameter. For this purpose, we write the coupling as u and we have from Eq. ([5]) 

F(v) = min[F(u) - u F'(u) + v F'{u)} . (12) 

u>0 

This is particularly useful in cases where f(s) is difficult to find explicitly. 

Another form of expression, useful for our present task, is obtained if we change the kinetic-energy parameter from s 
to r itself, by inverting the (monotone) function /(s) to define the associated X-function by 

K^(r)=s=(f- l of)(r). (13) 

Now the energy formula Eq. ([5]) becomes 

F(v) = min \K [f] {r) +v f(r)\ . (14) 

r>0 L 

A sleight of hand may be perceived here since K depends on /. However, we do now have a relation that has F on one 
side and / on the other: our goal is to invert this expression, to effect F — > /. We shall do this below by constructing a 
sequence of approximate if-functions which do not depend on /. 



B. Smooth transformations and envelope approximations 

In this section, we consider potential shapes f(r) that may be written as smooth transformations f(r) = g(h(r)) of a 
'basis potential' h(r). The idea is that we know the spectrum of — A + v h(r) and we try to exploit this to study the 
spectrum of —A + v f(r). When the transformation function g has definite convexity (g" does not change sign), the 
kinetic-potential formalism immediately allows us to derive energy bounds. This is a consequence of Jensen's inequality 
poj . which may be expressed in our context by the following: 

g is convex (g">0) (ip, g(h)tp) > g{{ip,hip)) , 

g is concave (g" < 0) (ip,g(h)ip) < g((ip,hi/})) . (15) 

More specifically, we have for the kinetic potentials 

g">0 =► f(s)>g(h(s)); g"<0 /» < g(h(s)) . (16) 

We can summarize these results by writing /(s) ~ g(h(s)) and remembering that the relation w indicates an inequality 
whenever g has definite convexity. The expression of these results in terms of X-functions is even simpler, for we have 

K [f] = t 1 o / « (g o h)- 1 o (g o h) = bT 1 o h = K [h] . (17) 

Thus, fa is the approximation we sought, that no longer depends on /. The corresponding energy bounds are 
provided by 

E = F(v) ~ min[s + vg(h(s)j^ = min 



(18) 
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C. The envelope inversion sequence 



We suppose that an eigenvalue E of H = — A+v f(r) is known as function E = F{v) of the coupling parameter v > 0. 
In some cases, such as the square well, the discrete eigenvalue may exist only for sufficiently large coupling, v > v\. The 
kinetic potential f(s) may be obtained by inverting the Legendre transformation in Eq. (|10|) . Thus 



F(v) = min[s + v /(s)] 



/(s) = max 



F(v) 



s>0 ' J v>vi 

We shall also need to invert the relation between F™ and K'- 71 ' by means of 

K{r) = ma,x[F(v) — v f(r)] ■ 



(19) 



(20) 



We begin with a seed potential shape (r) from which we generate a sequence (r)}^ =0 of improving potential 
approximations. The idea behind this sequence is that we search for a transformation g so that g(f^ (r)) is close to f(r) 
in the sense that the eigenvalue generated is close to F(v). The envelope approximation is used at each stage. The best 
transformation g^ at stage n is given by using the current potential approximation /^(r) as an envelope basis. We 
have: 



/ = 5 M o/ -N 



M — f o /W 



Thus 



j[n+l] _ [n] Q _ j Q ^[n] 



The resulting inversion algorithm may be summarized by the following: 
inversion algorithm 

/W( r ) fW (t) ) _ 



if [nI (r) = max F^(u) - u/ W (r) 



(r) = max 

t)>fi 



F(v) K^(r) 



(21) 
(22) 



The step /W (r) — >■ FW («) is effected by solving (- A + w /W) ip = Eip numerically for E = F^ (u). 



D. Uniqueness 

We consider now a singular potential f(r) of the form 

/(r) = i±± , where 5 (0) < , g'(r) > , (23) 
r 

and g(r) is not constant. Examples of this class of singular potential shapes f(r) are Yukawa g(r) = — e _or , Hulthen 
g(r) = — r/(e ar — 1), and linear-plus-Coulomb g(r) = —a+br 2 , with a, b > 0. With these assumptions, we have proved in 
Ref. [H the following 

Theorem 1 The potential shape f(r) in H = —A + vf(r) is uniquely determined by the ground-state energy function 
E = F(v). 



III. SPECTRAL DATA FROM MINKOWSKI-SPACE BETHE-SALPETER EQUATION 



A. The raw data 



In Table U we exhibit the binding energy E versus coupling v results from numerical solutions of the Bethe-Salpeter 
equation for a system of two scalar particles each of mass m, bound by single or multiple exchange of a scal ar p article of 
mass ii, computed by Carbonell and Karmanov in Refs. d, Table 1], d, Table 1], 7, Tables 1 and 2], and [l3[ Table 1]. 
For comparison, we add corresponding results of the Schrodinger equation with an interaction potential of Yukawa form 
Vir) — «exp(-/jr)/r, which, as shown in Appendix [XJ constitutes the nonrelativistic limit of the ladder BS quation. 
Interestingly, the nonrelativistic binding energies emerging from the Schrodinger equation with Yukawa potential seem 
to reproduce better the ladder-plus-cross-ladder approximation findings than the ones from mere ladder approximation. 
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TABLE I: Couplings v and binding energies E arising from Bethe-Salpeter equations in either ladder or ladder-plus-cross-ladder 
approximation for common mass m = 1 of the bound scalar bosons and mass /j = 0.15 or n = 0.5 of the exchanged scalar boson, 
or from the ladder-approximation nonrelativistic limit, the Schrodinger equation with Yukawa potential V(r) — —v exp(— /j r)/r. 



E 



Minkowski- Space Bethe-Salpeter Equation 



Schrodinger Equation 



Ladder 



Ladder + Cross-Ladder Yukawa Potential 



(j, = 0.15 % Z M M = °- 50 [Sri L3] M = 0.50 [6, 7} 



At = 0.50 



0.5716 

1.437 
2.100 
3.611 
5.315 



1.440 

2.01 

2.498 

3.251 

4.901 

6.712 



1.21 
1.62 
1.93 
2.42 
3.47 
4.56 



1.034 
1.285 
1.532 
1.848 
2.204 
2.918 



-0.01 
-0.05 
-0.10 
-0.20 
-0.50 
-1.00 



B. Exchange- mass dependence 

In this subsection, we demonstrate that the spectral data of Table Q] for fi — 0.5 can be obtained approximately from 
the corresponding data shown in Table[I]for /i = 0.15 by a scale change in the potential of a Schrodinger model. This is 
interesting if one expects to find that a potential such as the Yukawa V(r) — —v e~^ r /r would account approximately 
for the spectral dependence of the problem on the exchange mass fi. Let us consider a Schrodinger operator given by 

_A + U ZW _^ E = p, v s ^ (24) 
r 

where E is a discrete eigenvalue. A simple scaling argument applied to the operator 

H = -±A + vl^l (25) 
2m r 

shows that a corresponding discrete eigenvalue of H is given in terms of F(v) by the formula 

u 2 ( 2m v\ , . 

E =iL F — ■ (26) 



2m \ /i 

Thus, if we compare two different /i values, \i\ and fi2 with R = fj,i/fi2, and we write for the /xi case E\ = F\(v), then, 
under the scaling rule (|26| . for the //9 case we would have Ei = Fi(Rv)/R 2 . For our present problem, we have \i\ =0.5, 
jt*2 = 0.15, and therefore R = \ii = 10/3. Hence, given the second column of Table U expressed as E = E\ = F±(v), 
we would expect to generate data consistent with the first column of Table|T]by the formula E2 = (9/100) ^1(10^/3). 
Graphs of ^(w) and its approximation in terms of the scaled F±(v) are shown in Fig. [I] The scaling law seems to yield 
a rough approximation. 



IV. THE CONSTRUCTION OF EFFECTIVE POTENTIALS 



We now consider the BS spectral data collected in TableQ]and we use our inversion theory [3J] to answer the question 
what potential shape f(r) in the Schrodinger Hamiltonian H = — l/(2m) A + v f(r) would generate the corresponding 
binding energies E for the given values of the coupling parameter i>? For this purpose, we take m = mi m2j(rn\ + 7712) 
with mi = m,2 = 1- We adopt the inversion algorithm of Eq. (|22[) with the pure Coulomb seed potential J^^) = — 1/r. 
For the three sets of BS data in TableU our results are exhibited respectively in Figs. [1112] and SJ In each case, we first 
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FIG. 1: Approximate exchange-mass dependence from scaling arguments: Table[I]data ^(w) vs. scaled Table|I]data Fi(Rv)/R 2 . 

show a sequence of eight iterations and then depict the last potential iteration, (r), along with the eigenvalue curve 
F(v) of the Bethe-Salpeter data and the corresponding eigenvalue curve of the Hamiltonian H = — A + v /' 8 '(r). For 
comparison, Fig.[5]shows the reconstructed potentials for both ladder (L) and ladder-plus-cross-ladder (L + CL) kernels 
with an exchanged-particle mass \x — 0.5, along with the corresponding Yukawa-potential shape Y(r) = — exp(— /i r) / r. 
Judged by the eye, the Yukawa potential seems to be closer to the inversion output for the ladder-plus-cross-ladder case 
than for the mere ladder case. In view of the ordering of couplings v in Table HI this observation is no genuine surprise. 



V. FORM FACTORS 



In quantum mechanics, the three-dimensional form factor F(k) of a bound state described by its configuration-space 
wave function i/j(x) is nothing else but the Fourier transform of the corresponding charge density p(x) = \ip(x)\ 2 . For a 
given coupling strength v and common mass m of the two bound-state constituents, the reduced radial wave functions 
u(r) of s-states satisfy an ordinary differential equation that determines the associated binding-energy eigenvalues E: 

u"(r)+vf(r)u(r) =Eu(r) , u(0) = . (27) 

m 

It proves convenient to normalize the radial wave functions it(r) such that dr |u(r)| 2 = 1. In momentum space, any 
corresponding form factor F(k) is given by the Fourier-Bessel transform of the radial density u 2 (r), that is to say, by 

oo 

1 f sin(fcr) o, , „, . , . 

F(k) = — dr — i — '- u 2 ir) , F(0) = 1 . 28 
k J r 


Figures E] and [7] depict the ground-state form factors F(k) for v = 5 and m = 1 for each of the three potential shapes 
f(r) shown in Figs. [21 HI obtained by application of our geometrical inversion technique to the spectral data of Tablcfl] 
By comparing Figs.lBJa) and[7l we see that the quantum-mechanical form factors broaden with increasing mass fi of the 
exchange particle. Likewise, a comparison of the two plots in Fig. [6] reveals that the quantum-mechanical form factors 
broaden when taking into account higher-order corrections in the BS interaction kernel, that is to say, when working in 
the somewhat more sophisticated ladder-plus-cross-ladder approximation instead of in the naive ladder approximation. 
These form factors F(k) constitute the three-dimensional counterparts of the four-dimensional form factor F(Q 2 ) found 
from the ladder-plus-cross-ladder Minkowski-space BS equation with exchanged-boson mass fi = 0.5 in Refs. [8l ll(X[l3| . 
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o 



f(r) 




12 3 4 5 6 



(b) 



FIG. 2: Geometric inversion of spectral data of bound states of two scalar bosons described by a Minkowski-space Bethe-Salpeter 
equation in ladder approximation for mass fi = 0.5 of the exchanged particle [H-Qi EH • Graph (a) shows the first eight iterations 
f[n], n = 1, 2, . . . , 8, starting from the seed, while graph (b) shows the resulting potential shape f(r) and, in the inner graph, the 
corresponding spectral curve F(v) along with the input data. 
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ITr) 




(b) 



FIG. 3: Geometric inversion of spectral data of bound states of two scalar bosons described by a Minkowski-space Bethe-Salpeter 
equation in ladder-plus-cross-ladder approximation for mass fi = 0.5 of the exchanged particle [E01- Graph (a) shows the first 
eight iterations f[n], n = 1,2, ... ,8, starting from the seed, while graph (b) shows the resulting potential shape /(r) and, in the 
inner graph, the corresponding spectral curve F(v) along with the input data. 
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it.) 




(b) 



FIG. 4: Geometric inversion of spectral data of bound states of two scalar bosons described by a Minkowski-space Bethe-Salpeter 
equation in ladder approximation for mass /x = 0.15 of the exchanged particle [1,0, EH- Graph (a) shows the first eight iterations 
f[n], n = 1, 2, . . . , 8, starting from the seed, while graph (b) shows the resulting potential shape f(r) and, in the inner graph, the 
corresponding spectral curve F (v) along with the input data. 
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FIG. 5: Potential shapes /(r) from Figs.[2]and[3]resulting from geometric inversion of the ladder (L) and ladder-plus-cross-ladder 
(L + CL) Bethe-Salpeter findings [51471. 1 131] for the masses of two-scalar-boson bound states, with exchanged-boson mass /i = 0.5, 
compared with the corresponding Yukawa potential Y(r) = — exp(—fj,r)/r, arising in the nonrelativistic limit of the ladder case. 
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FIG. 6: Momentum-space form factor F(k) of the ground-state Schrodinger solution mimicking the results of the Bethe-Salpeter 
equation in ladder (a) and ladder-plus-cross-ladder (b) approximation with coupling v — 5 and exchanged-particle mass \i — 0.5. 
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1.2 



F(k) 




1 1 ' 1 ' ' ' ~ 1 ' 1 

2 4 6 8 10 

k 

FIG. 7: Momentum-space form factor F(k) of the ground-state Schrodinger solution mimicking the results of the Bethe-Salpeter 
equation in ladder approximation with coupling v = 5 and exchanged-particle mass fi = 0.15. 



VI. CONCLUSION 

For a quantum-mechanical description of relativistic systems in which a single particle is bound to a fixed center, one 
uses the Klein-Gordon equation or the Dirac equation. For the analysis of systems composed of more than one particle 
within quantum field theory, the Bethe-Salpeter formalism is required. In a series of papers, Carbonell et al. discussed 
the bound states of a system of two charged scalar bosons by means of the Bethe-Salpeter equation in Minkowski-space 
representation, in contrast to most studies of this kind which rely on the Euclidean-space formulation of this equation. 
These investigations report some numerical results for the binding energies, in a variety of cases, as functions E(v) of a 
coupling parameter v. What we have done in the present analysis is to employ a geometric spectral inversion theory to 
reconstruct, in each case, the potential shape f(r) in a Schrodinger model H = — A/ (2m) +v f(r) which would have the 
same energy curve E(v). As more complete quantum- field-theoretic spectral data becomes available, we shall be able to 
reveal more details of such spectrally-equivalent potential models. 
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Appendix A: Nonrelativistic Reduction of the Bethe-Salpeter Equation 



For the sake of completeness we briefly sketch how, by successive application of a sequence of simplifying assumptions 
and approximations, the Bethe-Salpeter equation for two bound-state constituents both of spin zero may be reduced to 
an equation of motion of Schrodinger form with all interactions represented by a static potential. Regarding kinematics, 
for a bound state of two particles discriminated by a label i = 1, 2 the relation between total momentum P and relative 
momentum p of the constituents, on the one hand, and the individual-particle momenta p\ , p-i , on the other hand, reads 



P = pi + p 2 



p = m pi - m P2 



Pi = Vi p + P 



P2 = m p - p 



where rjx, 772 denote two real parameters satisfying 771+772 = 1. In the center-of-momentum frame of some bound state of 
mass M = V^P 2 ", defined by P = P1+P2 = and therefore P = (M, 0), the individual-particle momenta pi, P2, become 

pi = r\ x M + p° , p\ = T] 2 M - p° , pi = p , p 2 = p ■ 

Our starting point of the nonrelativistic reduction is the Bethe-Salpeter equation in momentum-space representation 

®(p,P) = -^S l (p l ) J d i qK(p,q,P)<S>(q,P)S2(- P 2) ■ (Al) 

Its instantaneous approximation assumes that in the center-of-momentum frame of the bound state the Bethe-Salpeter 
interaction kernel K (p, q, P) depends only on the (initial and final) spatial relative momenta p, q: K(p,q,P) — K(p,q). 
Then, integrating over pq reduces the Bcthe Salpcter equation (|A1[) to a kind of instantaneous Bethe-Salpeter equation 



^(P) = 2^ / dPoSi(pi)S 2 (-p 2 ) 



d 3 q 
(2tt) 3 



K(p,q)4>{q) 



for the 'Salpeter amplitude' 4>{p) , defined as integral of P) over the time component po of the relative momentum p: 



(P) 



1 

2^ 



dp Q <£(p, P) 



Replacing the propagator Si(p) of bound-state constituent i by its free counterpart S^(p) entails the Salpeter equation 

d 3 q 



Hp) 



1 

2^ 



dpoS{ 0) ( Pl )si 0} (-p2) 



K(p,q)(t>(q) . 



(A2) 



The free Feynman propagator (pi) in momentum space of a scalar boson of mass mi and momentum pi is given by 



(0), 



p\ — m\ + i £ 



e|0 , 



1,2 



In terms of relativistic free-particle energies (p) = \Jp 2 + m? the unique 2 partial fraction decomposition of {pi ) is 



2Ei{p) 



p° - Ei(p) +ie p° + Ei(p) - is 



e|0 , 



1,2 



Evaluating the integral over the two propagators in Eq. (|A2I) by contour integration and Cauchy's residue theorem gives 



d Po s[ \ Pl )Sr(-P2) 



(0), 



2ni 



AE 1 {p)E 2 {p) 



1 



M - E x (p) - E 2 (p) M + E x (p) + E 2 (p) 



The three-dimensional reduction of the relativistically covariant Bethe-Salpeter equation (|Alj) to a Schrodinger-type 
equation is applicable to nearly nonrelativistic and weakly bound states composed of sufficiently heavy constituents. For 
such systems, the first contribution to the above integral over propagators may be assumed to dominate the second one, 



A different decomposition of the free scalar-boson propagator S± (pi) that is, however, not compatible with the Cauchy residue theorem is 



^ 0) ( Pi ) 



1 



1 



Pi 



Ei(p) + \e p1 + E t {p)-ie 
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1 1 

M - E^p) - E 2 {p) >> M + Ei(p) + E 2 (p) ' 
since in this situation M « E\ (p)+E 2 (p), such that the second term in the propagator integral may be safely neglected: 



4 Ei (p) E 2 (p) M - E 1 (p) - £ 2 (p) 
Adopting this standard approximation, one arrives at the reduced Salpeter equation for spin-0 bound-state constituents 

[E 1 (p)+E 3 (p)]i>(p)- —- \ f ■^K(p,q)4>{q)=M4>(p) . (A3) 

4E 1 (p)E 2 (p)J (2ny 

Moreover, assuming that, in the first term of the propagator integral, the factor [M—Ei(p)— E 2 (p)]~ 1 varies faster than 
the factor [E\ (p) E 2 (p)] _1 justifies the substitution Ei(p) m^ in the denominator of the interaction term in Eq. (|A3|) . 
Finally, the nonrelativistic expansion Ei(p) w m^+p 2 / (2 m^) of the free energies leads to the Schrodinger- type equation 

mi +m 2 + ^— + HP) - I K & 9) #9) = M • ( A4 ) 

2 mi lm 2 ) 4mim 2 J (Zir)° 

The Schrodinger equation governing the dynamics of two particles, of masses mi , 1112, interacting via a potential V(a?) 
involving their relative coordinate x = Xi—x 2 , in a bound state of mass M reads, in configuration-space representation, 



mi + m 2 - h V(x) 

2 mi 2m 2 



ip(x) = Mip(x) 



Upon introduction, for configuration-space wave function ip(x) and interaction potential V(x), their Fourier transforms 

the latter proving to be a very convenient choice, the Schrodinger equation becomes in momentum-space representation 

mi+m 2 + + #(p) + / y(p - q ) = M ^(p) . (A5) 

2 mi 2m 2 / J ( 2 7r) 

Assuming the kernel K(p, q) to be of convolution type, i.e., K(p, q) — K(p—q), the comparison of the reduced Salpeter 
equation in nonrelativistic limit (|A4[) with the momentum-space Schrodinger equation (|A5I) allows for the identification 



V( p -q) = -—± K{p-q) «=> F(p)=_^ Jjf(p) . 

Amim 2 4mim 2 

As an illustration of this relationship, let us demonstrate how the Yukawa potential arises from single-boson exchange 
between our two bound-state constituents i — 1,2. Let <?i denote the interaction strength (having the mass dimension 1) 
of the three-boson coupling of the spin-0 bound-state constituent i (of mass mj to some spin-0 force mediator of mass [i. 
Introducing the momentum transfer k = p—q, the resulting one-boson exchange contribution to the interaction kernel is 

i (i^i) (152) 



iK(p,q,P)=iK(k) 

r approx; 
K(k) = 



k 2 — (j? 

contenting oneself with this form entails the ladder approximation to Eq. (|Al|) . In instantaneous limit, this kernel reads 

5i52 



2 



k 2 + n 

As consequence of the spherical symmetry of K (fe), the Fourier transformation of V(k) yields the spherically symmetric 
configuration-space Yukawa potential V(x) = V(r), r = \x\. For convenience, we represent it in the form V(r) = v f(r): 

V ( x ) = / T^ZvT e V ( k > = - / ToZvT , =- ,a — = V(r) = vf(r), f(r) 



(2tt) 3 J (27r) 3 4mim 2 167rmiTO2 r r 

Hence our coupling constant v is related to the mass and coupling parameters of the underlying quantum field theory by 

9i 9i 
v = . 

167rmi m 2 

Especially, for identical bound-state constituents, clearly satisfying mi = m 2 — m and gi — g 2 = g, this result becomes 

9 2 
167r m 2 
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